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Abstract. A consortium of universities has recently been formed with the goal of using the decommissioned telecom- 
O . munications infrastructure at the Goonhilly Earth Station in Cornwall, UK, for astronomical purposes. One partic- 

ular goal is the introduction of one or more of the ~30-metre parabolic antennas into the existing e-MERLIN radio 



CZ2 ' interferometer. This article introduces this scheme and presents some simulations which quantify the improvements 

^ that would be brought to the e-MERLIN system. These include an approximate doubling of the spatial resolution 

of the array, an increase in its N-S extent with strong implications for imaging the most well-studied equatorial 
fields, accessible to ESO facilities including ALMA. It also increases the overlap between the e-MERLIN array and 
^ ■ the European VLBI Network. We also discuss briefiy some niche science areas in which an e-MERLIN array which 

\l ' included a receptor at Goonhilly would be potentially world-leading, in addition to enhancing the existing potential 

, of e-MERLIN in its role as a Square Kilometer Array pathfinder instrument. 

1. Introduction In this short article we present some simulations to 

rr., ^ , .„ , r, ■ / ^ ^ , ^^^^^-^t , demonstiate thc Improved peiformance of au e-MERLIN 
The Goonhilly Earth Station (GES; lat = 50.0504 N, Ion = u- u • i j \. \ i^r i j- 

^ . , , , / N , ,. which mcludes an antenna at Goonhilly. We also discuss 

5.1835 W) is home to three large (~30-metre) parabolic . ..n i- i- i u u -j i r i, 

' , . , • . ,1 1 ., r 1 • some scientific applications which would be ideal tor such 

antennas which were oriemally built tor telecommumca- . , , rpi • i- i • i r u ^ 

° ■' an instrument, ihis article is also a follow- up to a previ- 
tions. GES was in fact the world s first satellite commu- i,- u j- u cj. ^ iu • i • ^ 

• ^ ous paper which discusses the benefits ot the inclusion or 

k> ■ nication station, its original parabolic antenna being used ii oc i r a a r^u-iu li r^u j. 

, . , ' , . .... , . the 25-metre parabolic antenna at Ghilboiton Observatory 

, to receive the farst transatlantic television broadcast via j j. i onno^ 

i3 ' ^ , ,,. . , , (Heywood et al., 2008). 

■ the ielstar satellite approximately hali a century ago. 

Since most satellite communication and international 
telecommunication needs are now met using a combina- . The e-MERLIN array which is described in more de- 
tion of much smaller dishes and undersea cables, several ^^^^ Section 2, serves as the UK's Square Kilometre 
of which terminate at the Goonhilly site, the satellite op- ^^^^ay (SKA) pathfinder facility It has a unique science 
erations using the large dishes at Goonhilly were ceased scope due to its "intermediate" baseline lengths (i.e. oc- 
in 2008. The Consortium of Universities for Goonhilly cupying the gap between those of the EVLA and VLBI 
Astronomy (CUGA0; consisting of the universities of networks such as the VLB A and EVN). Its scientific ex- 
Exeter, Hertfordshire, Leeds, Manchester and Oxford) has ploitation can thus potentially strongly infiuence the even- 
recently formed in order to use the existing infrastructure ^ual placement of the dishes in Phase-I of the SKA (SKAi ; 
at Goonhilly for astronomical purposes. One such applica- Garrett et al., 2010). Goonhilly may be particularly im- 
tion is the connection of one or more of the large parabolic Portant in this respect as it enhances the resolution of e- 
antennas into the e-MERLIN array MERLIN enough to allow the exploration of specific SKA 
Key Science aspects which are certain to require such in- 



^ [http : //www ■ ast . leeds . ac . iik/~mgh/CUGA . htm | termediate baseline lengths, but for which there is little 
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Fig. 1. The locations of the existing telescopes of the e- 
MERLIN array and the location of the Goonhilly Earth 
Station. The additional long baselines provided by a sta- 
tion at Goonhilly are immediately apparent. 

direct observational evidence concerning just what range 
of baselines is needed. 

The high cost of data transport from outlying dishes 
(i.e. the receptors which form the intermediate and longest 
baselines in the SKA) means that e-MERLIN will also 
inform on the trade-off between science return and project 
cost for SKAi. 

2. e-MERLIN 

The e-MERLIlsd array, operated by the University of 
Manchester at Jodrell Bank Observatory, is a radio in- 
terferometer consisting of seven radio telescopes situated 
around the UK, with a current maximum baseline of ap- 
proximately 217 km. Upgrades to the original MERLIN 
array include new receivers and a high-bandwidth fibre 
link to transport the data back to a new digital correla- 
tor. 

The upgraded system allows e-MERLIN to deliver 
micro- Jansky sensitivity with up to 4 GHz of instanta- 
neous bandwidth at L- (1.3 - 1.8 GHz), C- (4 - 8 GHz) 
and K-band (22 - 24 GHz), with an angular resolution of 
10 - 150 milliarcseconds. 

The locations of the existing e-MERLIN telescopes are 
shown by the white markers on Figure [U and also marked 
on this map in yellow are the locations of the GES and 
Chilbolton Observatory. Note that Jodrell Bank hosts two 
antennas: the 76-metre Lovell Telescope and the 25-metre 
Mark-II. The extra long baselines formed by connecting a 
receptor at Goonhilly to each of the existing e-MERLIN 
telescopes should be immediately apparent. 

^ http : //www. e-merlin. ac .uk) 



The surface accuracy of the Lovell telescope precludes 
its inclusion in K-band observations, and similarly the 
telescope at Defford cannot be used because of the spac- 
ings of its mesh refiector. 

Note that even if both the Lovell and Mark-II tele- 
scopes are included in an e-MERLIN observation there is 
still one spare correlator input which could be fed by an 
additional antenna. 

3. Extra baselines and improved uii-coverage 

The GES has three large antennas on site: Goonhilly-1 (or 
"Arthur"; 25.9 m diameter), Goonhilly-3 (or "Guinevere"; 
29.6 m diameter) and Goonhilly-6 (or "Merlin", no rela- 
tion; 32 m diameter). We assume that the proposed ex- 
tension of e-MERLIN would initially involve the use of 
Goonhilly-1 at L- and C-bands, and the eventual use of 
Goonhilly-3 at C- and K-bands. 

In an interferometric observation each baseline as pro- 
jected on the sky measures a single Fourier component of 
the sky brightness distribution at a particular observing 
frequency. A single correlation product between two an- 
tennas is the atomic unit of an interferometric observation, 
a "visibility" , and this complex number and its conjugate 
occupy distinct points in the uw-plane. The longest base- 
line of an array occupies the outermost region of the uv- 
plane and provides the highest resolution. Similarly the 
shortest baseline samples the largest spatial scales and 
occupies the inner-most region. The shortest baseline also 
determines the spatial scale at which objects begin to get 
"resolved out" . 

In addition to increasing the sensitivity via an increase 
in collecting area, adding more antennas to an array in- 
creases the range of spatial scales it can sample by virtue 
of the the extra baselines which are created, thus adding 
more points to the ww-plane coverage. Earth-rotation in- 
creases the uu-coverage further by changing the projection 
of the baselines on the sky with time, causing the baselines 
to trace arcs on the uv plot. Multi- frequency synthesis 
can also be used to increase the spread of measurements 
in the uw-plane, by observing over a broad bandwidth, 
which manifests itself as a radial spread in the uv plot. 

We can demonstrate the improved uw-plane coverage 
that the introduction of Goonhilly would bring by simu- 
lating an e-MERLIN observation with an additional an- 
tenna at the location of Goonhilly-1. This is accomplished 
by generating a set of simulated visibilities^ The antenna 
table is constructed using the locations of antennas from 
an actual e-MERLIN observation with an additional entry 
for the Goonhilly-1 antenna. 

We simulate a full earth-rotation L-band observation 
with 256 channels between 1.2 and 1.712 GHz for three 
representative declinations, and an antenna elevation cut- 
off of 8 degrees is also imposed. 



Simulations are performed using the sm tool within the 
CAS A ( http : //casa ■ nrao . edu J package. 
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Fig. 2. Four full synthesis uu-coverage plots for the e-MERLIN array with a frequency range of 1.2 to 1.712 GHz, for 
four source declinations as indicated on each panel. The extra baselines provided by the Goonhilly station are shown 
in red. 



Figure[2]shows the ww-plane coverage of the e-MERLIN 
array with the inclusion of the Goonhilly antenna. This 
figure is constructed by plotting the resulting u and v co- 
ordinates from the simulated visibilities at four represen- 
tative declinations and colouring in red any baseline that 
is formed with the Goonhilly- 1 antenna. 

These plots clearly demonstrate the increased and im- 
proved spatial frequency sampling and resolution advan- 
tages of an extended e-MERLIN array. Of particular note 
is the upper panel showing the declination = case. When 
observing the celestial equator with an interferometer the 
baseline tracks are in a purely E-W direction, generally 
resulting in undesirable N-S structure in the dirty beam 
(or point-spread function). Any N-S extension of the array 
itself goes some way to alleviating this problem. 

Note that for the simulations presented here we in- 
clude the Mark-II but not the Lovell Telescope. Generally 
only one of these two telescopes is included in the interfer- 
ometer, and inclusion of the Lovell instead of the Mark-II 



would make little difference to the w?;-plane coverage as 
the telescopes are located very close to one another. The 
main effect the Lovell Telescope has when included in e- 
MERLIN (at L- and C-band) is to form a series of very 
sensitive baselines due to its 76-metre diameter, increasing 
the sensitivity of the array by a factor of 2-3. 

Although the inclusion of the Lovell has implications 
for wide-field observing due to its comparatively small 
field-of-view, the proposed e-MERLIN Legacy Surveys will 
generally benefit from its inclusion, and such programmes 
especially those targeting faint objects at high resolution, 
would be enhanced greatly by the very sensitive long- 
baseline formed by Goonhilly and the Lovell (Klockner 
et al., 2011, in prep.) 

4. Improved beam shapes 

The visibilities used to generate the uv plots in the pre- 
vious section are imaged with natural weighting to deter- 
mine the dirty beam (i.e. point-spread function) of a given 
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Table 1. Major axis, minor axis and position angle of 
the Gaussians fitted to the central lobes of the naturally- 
weighted dirty beams for full-track observations at four 
declinations at L-band. These effectively correspond to the 
resolution of the array. Also included in the table are the 
ratios of the major to minor axis for each scenario. The 
closer this number is to unity the more circular the beam, 
and Goonhilly provides vastly more circular beams for low 
— Dec — observations. (eM — the existing e-MERLIN ar- 
ray; eM-l-G includes Goonhilly-1). 



Declination 






PA 


Bmaj / 


(and array) 


(arcsec) 


(arcsec) 


(deg) 


Bruin 


5 = -20° (eM) 


0.688 


0.159 


12.72 


4.33 


5 = -20° (eM+G) 


0.352 


0.118 


163.482 


2.98 


<5 = 0° (eM) 


0.393 


0.166 


22.28 


2.37 


5 = 0° (eM+G) 


0.195 


0.116 


-29.65 


1.68 


5 = -f 30° (eM) 


0.225 


0.186 


23.45 


1.21 


(5 = -H30° (eM-hG) 


0.135 


0.120 


134.39 


1.13 


5 = -f60° (eM) 


0.200 


0.183 


0.74 


1.09 


5 = +60° (eM+G) 


0.133 


0.126 


-1.06 


1.06 



observatior0 The generation of these beams naturally in- 
corporates multi-frequency synthesis. 

A two-dimensional Gaussian is fitted to the central 
lobe of each of the generated beam patterns and the results 
of this procedure are presented in Table [T] In addition to 
the major and minor axes of the Gaussian their ratio is 
also tabulated. This provides a measure of the "circular- 
ity" of the beam. 

Having a point-spread function which is more circular 
for equatorial fields will benefit survey science in general, 
as many of the best-studied multiwavelength fields are sit- 
uated on or near the celestial equator, and it is likely to be 
the site of the deepest future ALMA observations. Circular 
beams are also of particular benefit to weak lensing exper- 
iments. In order to overcome the systematic limits in deep, 
wide-area weak lensing surveys, and make the first steps 
towards exploiting weak lensing with the SKA (Blake et 
al., 2004), there is an urgent need to extend existing first 
attempts at radio-optical weak lensing surveys (Patel et 
al., 2010) to cover larger-area weak-lensing surveys, typi- 
cally situated at or below the celestial equator. 



* When imaging an interferometric observation the visibili- 
ties are gridded and then Fourier transformed. The Mw-plane 
coverage and the dirty beam are directly coupled, in that 
they are Fourier transforms of one-another. Natural weight- 
ing means the gridded visibilities are weighted according to 
the density of points in a uv cell. Since there are always more 
short baselines than long ones naturally-weighted visibihties 
result in a lower resolution than that returned by different 
weighting functions. Uniform weighting apphes equal weight 
to all visibilities and generally results in the highest resolution 
at the expense of imaging sensitivity. 



As mentioned previously, the baselines during a dec- 
lination = observation move through the ww-plane in 
a purely E-W direction as can be seen in Figure [2J This 
generally results in unfavourable N-S sidelobes in the dirty 
beam, and increasing the N-S extent of e-MERLIN by in- 
cluding Goonhilly mitigates this problem significantly. 



5. Applications 

Here we briefly discuss some e-MERLIN science cases for 
which the enhanced performance of the instrument would 
be particularly valuable. 



5.1. Strong gravitational lenses (primarily L-band) 

The concomitant advent of an efficient means of find- 
ing strongly-gravitationally-lensed systems (HERSCHEL; 



|Negrello et al., 2010 1 and a means of following these up 
at high spatial resolution (e-MERLIN; Deane et al. 2011, 
in prep.) opens up some fascinating new possibilities. 
Goonhilly provides e-MERLIN with "HST-like" resolution 
even at L-band where the lensed radio sources are typically 



brightest. As emphasised by Dalai & Kochanek (2002) 
and Vegetti fc Koopmans (2009)[ such resolution is suf- 
ficient and necessary to adequately probe the dark mat- 
ter distribution of quad-system gravitational lenses down 
to sub-galactic scales. Critically, the new HERSCHEL / 
e-MERLIN method will target lensed sources that are 
intrinsically extended enough to avoid microlensing ef- 
fects, but compact enough to be noticeably perturbed by 
small-scale substructure in the dark matter. By measur- 
ing the mass fraction in sub-structure in a large sample of 
new gravitational lenses, e-MERLIN with Goonhilly could 
prove to be the current-generation facility best-suited to 
testing the predictions of the CDM model on mass scales 
well below that of the Milky Way satellites. 

5.2. European VLBI (primarily L- and C-band) 

Forming baselines between the existing e-MERLIN tele- 
scopes and Goonhilly will also improve the overlap be- 
tween e-MERLIN and the European VLBI Network, which 
has a minimum baseline length of ~270 km (Effelsberg 
- WesterborlQ. Combined EVN+MERLIN observations 
are currently routine, where the shorter baselines of 
MERLIN are required to resolve larger scale structure. 
Adding a Goonhilly receptor to the e-MERLIN array will 
introduce overlaps in baseline lengths between these two 
arrays. Goonhilly would also increase the western extent 
of the EVN. For a more detailed discussion of this aspect 
see the article by Klockner et al. (2011, in prep.). 



^ Note that 12 of the 14 elements which form the Westerbork 
telescope are being fitted with the APERTIF focal-plane array 
system in 2012. Since only two dishes will remain available for 
VLBI work above L-band, the EVN will have greatly reduced 
C-band sensitivity on this baseline. 
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Fig. 3. Simulations of observations of a CO disk at a redshift of 4.023. Simply to illustrate the channelised source 
morphology, the left-hand column shows thumbnail images of a selection of numbered channels from the input model, 
which is the molecular gas disk of a simulated galaxy at z = 4.023 as described by the S'^-SAX simulation. The four, 
large colour frames show the frequency-averaged cube emission for the input model (upper- left), an ALMA Band 4 
observation of CO(7-6) (upper-right), a current e-MERLIN observation of CO(l-O) at K-band (lower-right) and a 
second e-MERLIN observation but with the Goonhilly antenna included (lower- left). Each simulated observation has 
24-hours of on-source time, and colour scales run linearly from the minimum to the maximum pixel value. Contour 
levels are (-4, -2, 2, 4, 8, 16, 32...) x 0.01 mJy per beam for the three simulated images and (-4, -2, 2, 4, 8, 16, 32...) x 
0.01 mJy per pixel for the input model. This figure shows that the resolution provided by an e-MERLIN array which 
includes a Goonhilly antenna is sufficient and necessary to resolve and recover the disk orientation in a typical z^4 
massive galaxy. 



5.3. Massive star formation (primarily C- and K-band) 

To understand the complex dynamics driven by gravity, 
rotation, radiation and magnetic fields during the process 
of massive star formation high resolution probes of both 
the infalling molecular material and outflowing ionized gas 
are needed. ALMA will provide the required sensitivity to 
molecular gas at a resolution of around 10-100 mas which 
corresponds to 30-300 AU at a typical distance of 3 kpc 
for the nearest samples. Only e-MERLIN can provide the 
matching resolution and sensitivity at cm wavelengths to 
probe the ionized jets and winds that are an integral part 
of the star formation process. 

Significant progress will be made only through joint 
high resolution studies of a common sample which neces- 



sitates e-MERLIN with Goonhilly being able to observe 
equatorial objects also accessible to ALMA. Furthermore, 
most massive star formation takes place in the in- 
ner Galaxy (e.g. Urquhart et al., 2008) and repre- 
sentative samples require observations in the range - 
20° < Dec < -1-2°. Most of the best examples of jets 
from massive protostars are found in this region (e.g. 
Carrasco-Gonzalez et al., 2010; Gibb et al., 2003) and only 
e-MERLIN with Goonhilly can provide the factor of ten 
increase in resolution over current studies to probe the 
inner regions where the flows are being driven and coUi- 
mated. 
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5.4. Resolving molecular gas disks at high redshift 
(primarily K-band) 

The ^^CO molecule is the best proxy for molecular hy- 
drogen and its J(l-O) transition {lyrest = 115.27 GHz) is 
redshifted in to the K-band of e-MERLIN at z~4. High- 
redshift detections of molecular gas typically target sys- 
tems which are undergoing vigourous episodes of AGN ac- 
tivity (e.g. quasars) or star-formation (e.g. sub-mm galax- 
ies). ALMA will be the natural choice of instrument to 
target the higher-order transitions in such galaxies how- 
ever it is incapable of observing the lowest-order transition 
which provides the most robust estimate of the total mass 
of molecular gas. 

The S^-SAX simulation ( [Obreschkow et al, (2009)] 
and references therein) makes strong predictions about 
the cosmological evolution of molecular gas disks. We can 
demonstrate the advantage of adding Goonhilly to this ar- 
ray by selecting a massive object at z = 4.023 from the 
S^-SAX databas^l and simulating a line observation on 
the celestial equator. Galaxies are extracted from the sim- 
ulation catalogue within a random volume of space defined 
by the field-of-view of a 25-metre antenna along the spatial 
axes and the maximum K-band bandwidth along the red- 
shift axis. A model sky cube is then generated (Levrier et 
al., 2009) for the most massive galaxy within this volume 
showing the CO (1-0) emission from its molecular disk, red- 
shifted to 22.948 GHz, covered by 32 channels of 21.4 MHz 
each. For comparison we also simulate an ALMA observa- 
tion of the CO(7-6) emission {I'rest = 806.89 GHz) which 
is redshifted into Band 4 of ALMA at 160.63 GHz. A se- 
lection of numbered channels from the CO (1-0) cube are 
shown in the left-hand column of Figure |3] to illustrate the 
morphology of the source. 

These two model datacubes are then inverted to gener- 
ate three sets of visibilities, namely for an e-MERLIN ob- 
servation both with and without Goonhilly, and an ALMA 
observation using a current proposed layout for its most 
extended configuration. On-source time in all cases is 24 
hours, and noise is added accordingly. The simulations 
are performed in spectral-line mode and imaged with uni- 
form weighting for maximum resolution. The final, decon- 
volved, frequency-averaged images are presented in Figure 
[3l Clockwise from the upper-left these are the input model, 
the ALMA simulation, the e-MERLIN simulation and the 
e-MERLIN simulation including Goonhilly. Please refer to 
the caption for contour and colour-scale details. 

One of the fundamental requirements for ALMA is ex- 
cellent image fidelity and ALMA faithfully recovers the 
disk structure in our simulation. The current six-element 
e-MERLIN array barely resolves the source and because 
of the beamshape it mis-informs as to the disk alignment. 
The image produced by e-MERLIN including Goonhilly 
comfortably resolves the gas disk. For the three simu- 
lated cases this would be the most valuable observation 
in terms of constraining the gas kinematics by spatially 
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and spectrally resolving this lowest CO transition in this 
high-redshift object. Note that the resolution of an ex- 
panded e-MERLIN array at K-band would exceed that of 
ALMA for any frequency below Band 7. 

Even for galaxies which are barely resolved (e.g. >3 
independent points) it is possible to derive a crude rota- 
tion curve from the red- and blue-shifted measurements 
either side of a mean. Such observations can be combined 
with models for gas kinematics in clumpy disks to de- 
rive estimates of gas and dark-matter masses in distant 
z~1.5 galaxies (e.g. Daddi et al., 2010). An expanded e- 
MERLIN array may well be the best facility for extend- 
ing such studies to z~4 since obtaining such high spatial 
resolution with ALMA would require targeting of CO at 
impractically high-order J-transitions. 

Note also that high-redshift CO is another area 
in which an extended e-MERLIN would be a unique 
pathfinder instrument for influencing the design of the 
SKA. New and upgraded cm- wave facilities such as 
the EVLA (e.g. Riechers, 2010), Greenbank Telescope 
(e.g. Harris, 2010) and eventually MeerKAT (Heywood 
et al., 2011, these proceedings) are becoming increasingly 
important for studies of the low-order transitions of CO 
at high-z. 

A realisation of the SKA with e-MERLIN-like base- 
lines will be capable not only of detecting and resolving 
such galaxies, but doing so with extremely high detection 
rates via the survey speed offered by the extreme sensitiv- 
ity and wider field of view (due to the smaller dish sizes). 

6. Conclusions 

Introducing one of the 30-metre class antennas at the 
Goonhilly Earth Station into e-MERLIN would double the 
already superb resolving power of the array by increasing 
the maximum baseline from 217 km to 441 km. 

During the course of an observation an interferometer 
samples a range of spatial scales governed by the baselines 
between antennas as projected onto a sky plane perpen- 
dicular to the observing direction (thus defining the uv- 
coverage), and including Goonhilly results in significant 
improvements in this aspect. Of particular note is how 
the increased N-S extent of the array improves its per- 
formance during low-declination observations. The point- 
spread function of an aperture synthesis observation is es- 
sentially the Fourier transform of the uw-coverage, and the 
improved ww-coverage when Goonhilly is included results 
in correspondingly lower PSF sidelobes. 

An antenna at Goonhilly also increases the overlap in 
baseline lengths between e-MERLIN and the EVN, and 
would also increase the E-W extent of the latter (Klockner 
et al., 2011, in prep.). 

We have highlighted some science cases where a 
Goonhilly-enhanced e-MERLIN would be an extremely 
powerful instrument, including resolving molecular gas 
at high-redshift with resolution exceeding that of all 
but the highest ALMA bands, and follow-up observa- 
tions of newly-discovered strong gravitational lens sys- 
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terns. The latter is particularly noteworthy as an expanded 
e-MERLIN may well be the best contemporary facility for 
testing the predictions of the CDM model on small mass- 
scales. 
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